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Introduction 


Fractional calculus and its applications have attracted the attention of academics in recent years, not only in 
mathematics but also in a wide range of other sciences, including physics, thermodynamics, engineering, 
economics, etc. Electrical, electrochemistry, statistics, and probability are just a few of the fields where 
fractional calculus is used extensively. Furthermore, many cosmic events that conventional differential 
equations are unable to describe can be described by fractional differential equations [1-3]. 


One of the key techniques for estimating the solution to differential equations is the homotopy permutation 
approach, which has been the subject of several studies. One of the integrals transforms utilized in the 
solution of differential equations is the natural transform. To discover the approximate solution to differential 
equations, we may now combine the homotopy permutation with the natural transform. [4-6]. 


Researchers in this discipline have thrived from their initial steps, due to numerous applications in physics, 
economics, engineering, and biological sciences. Several complicated and efficient techniques, such as [7- 
34], have been created and developed to solve fractional differential equations. The goal of this paper is to 
provide both approaches for solving FPDEs using the fractional Atangana-Baleanu-Caputo operator 
(FABCO), NVIM, and NHPM. 
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In this paper, we apply NVIM, and NHPM to find solution of fractional differential equations with the 
fractional operator Atangana-Baleanu. The order of the paper is as follows: The basic definitions for calculus 
and fractional integration are presented in section 2, the methods used are analyzed in section 3, many 
examples are given that explain the effectiveness of the method proposed in section 4, and finally, the 
conclusion is provided in section 5. 


Basic concepts 
Definition 1. [35] The Atangana-Baleanu-Caputo operator when f € H*(e1, €2),€1 > €2 and forO << S <1 


1S 


T oy S 
ABCDS F(t) == fF fx) Es (-e-2)8)dx, (1) 
where B(S) is a function such that 8(0) = B(1) = 1. 


Definition 2. [36] Assume set of function, 


Il 
A= fig AM, a1,a2 > 0,|f(t)| < Me“,t € (-1)/ x [0, a} 
Then the natural transform over A set 
N(FC@E)) = Rus) = fe fut) de, (2) 
0 
The natural transform yields the Laplace transform via the following formula[37], 
lax Los 
aoe —st/u <a = 
R(u,s) == | est!" F(t) dt = —F (-). (3) 
0 
From [38] and Eq.(3) , we get this relation 
B(S) 1 
N (4°°D$ (2) = ——_, (x (us) - =£(0)) (4) 
1-s+5(=) 
Definition 3. [37] The inverse natural transform of a function is defined by 
p+ico 
1 st 
N-1(R(u,s)) = f(t) = a i eu R(u,s)dt, (5) 
p—ico 


where s and u are natural transform variables and p is a real constant. 
Analysis of NVIM 

Let’s considering FPDEs with FABCO 

AEDS O(u,t) + L(O(7)) +M(OU7)) =f), (6) 


with initial condition ®(u,0) = ®y(u) , where “8°DS$ is the FABCO, L is the linear operator, M is the 
nonlinear operator and f (uu, T) is a source term. 


We obtain by applying the natural transform to Eq.(6) with the stated initial condition 


(co) -20(u,0)) = NIFH) -L@)-M@], 


B(S) 
imsus(Zy 
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by substituting initial condition of Eq.(8), 


® =10,(4) 2) wip -L@)-M@), 


B(S) 


appling variation intretion method 


1-s+s(¥) 


= = = 1 
Ongi = Oy, +A Py — = Pow) + B(S) 


NIL(®,) +M(®,)- fo] ]}, (©) 


Note that 4 is the Lagrange multiplier, and since 0 < S < 1, hence A = —1, we obtain after applying the 
inverse of the natural transform to both sides of the equation 


_,{t-5+8(§) 
Ory, = Po(u) —N sis) Nn) + M(®,) —f(u,t)] }, (10) 


where we have the initial iteration is ®(u, 0) = ®y(u), and hence 
D(u,T) = limy soo Px (UT). 
Analysis of NHPM 
Subject to the provided initial condition, apply the natural transform to Eq.(6). 
28, (w(u,s) - 200) = NGO) -Lle@AI-Mlo@OD, AD 


1-5+5(“) 


by substituting initial condition of Eq.(11), 


1-s+s(*)’ 


~~ NM (L[v] + M [vl - 9), (12) 


# = v(x) - 
re) 
Taking the inverse of the natural transform and applying it to both sides of the equation (12), 


isis) 


1 s+s(%)° 
v=) (x) + N7-1| —~“+ N(g) | -N7 a6) 


B(S) N(Llv]+M [v]) }], (43) 


applying homotopy permutation method, 


U(4,t) = Yn=0P"Un(% 4), NU, 4)] = Xn=0P"Hn(w), = (14) 


where 


Hy (04) U2, V3) 0 Un) = i(x,t))]p-0,n = 0,1,2,. 


n! sah é 
Substituting Eq.(14) into Eq.(13) gives _ the result that, 


1-545(2) 


SN ro P" LO] + Lr-oP"Hn(w)) | |, (15) 


Yn=0 2? Un (x,t) = G(x,t) —@p Nn 


on comparing both sides of Eq.(15), the following result is obtained, 
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Pp: (x,t) = G(%,t), 


i=s25(4) 


pi:0,(%,t) =-N BS) 


N (Llwo] + Ho(w)) 


1-s+s(%) 


pv, («,t) =—-N1 BS) 


N(Llon-1] + Hy-1(v)) , (16) 


using the parameter y, we expand the solution in the following form 


(oe) 


u(x,t) = > pv, (4%), (17) 


n=0 


choose p = | results in the solution of Eq.(11) 


u(x,t) = lim Ds pv, (4, Ft) = » V(x, €). (18) 
PP 
n=0 n=0 


Application 


In this part, two nonlinear equations and a linear system will be solved using the NVIM and NHPM 
techniques, with the assumption that B(S) = 1. 


Example 1. Consider the nonlinear equation with the FABCO 

ABCD S (yt) = -< (=e = u) ® = Beet) 
Ou\u Op? 

depending on the initial condition ®(u, 0) = pu? 

Below we present the NVIM, 

Applying the NVIM to Eq.(19), we get 


or us\ 12 12 
Ong = -N[ (1-5+5(2) N (= Onn — 5 OR ~ (@2 Yun + Pn) ), (20) 


now, we find the approximate solutions as, 


uy> : 
oan (1-545) )waa)) = ‘(2 D+sgop) 
uy> 
®,=y?+N 1 (1-5+5(Z) )n (o- +85) 


a 2S 
-#(1+a-s@- StGias) -so= Nesxyt Sap) (21) 
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and so on. 


Therefore, the series solution ®(, T) of Eq.(19) is given by 


a 


25 
(yu, 2) = p? ((S? — 35 + 3) + (35 — 257) + S*@ ts), (22) 


when choosing S = 1 in Eq.(22), it becomes 


Saree (i4r4t4) (23) 


Now, we present the NHPM, 
applying the NHPM to Eq.(19), we get 


din=0 p Ay _ 2n=0 p' Bn 
+ din=07 Cn + Yn=0P" On 
by comparing both sides of the Eq.(24), the following result is obtained, 


Tp"), = p2— pn (2 = as (<)) | I) (24) 


D:D = Lu’, 


S 
p:®,=N7! (2 -s+5(*) ) Mp? Ay ~ p By, + pCy + pol} 


J 


AL Ss 
p20, = N-! (: -s+8(—) Jarier, — p'B, + pC, + 9'0)] 


by the above algorithms, 
Dy =H, 

2 a 
0,=wu (1-stsa) 

1 125 
i= (« — 25 + 87) + (25 - sant 8? aa) 
and so on. 

Therefore, the series solution®(y, T) of Eq. (19) is given by 


S T2725 


O(u,T) = we? cs — 38 + 87) + (3S - 28) apt Srasapt | (25) 


If we put S > 1 in Eq.(25) , we get the approximate and exact solution 


a. 
oud =w (14 t+ oe), (26) 


Ultimately, the exact solution of Eq.(19), 
OCT) =pre*.. (27) 
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Example 2. Consider the nonlinear equation with the FABCO 
1 
Dz ®(yu,t) +5(®"),-O+ O° =0,0<6 <1 (28) 


with the initial condition ®(u,0) = e". 
Below we present the NVIM, 


using algorithm of the method, we get 


ny, =e — NU (( Sa ()) N (o + 5(®), — »)) (29) 


Now, we discover the approximate solutions as follows:®) = e+, 


S 
oer -((1-545()) N(- oH)=e (2- (2-stsz5 55) 
S 
®, =e — (1-545) ‘\w we t ~_ 
T2s 


thus, the approximate solution of Eq.(28) can be written, 


a 25 
oun = (0 3S + $7) + (3S - 2S apt ‘apt ) (31) 


when choosing S = 1 in Eq.(31), it becomes 
a 


O(u,T) =e it ‘), (32) 
Now, we present the NHPM, 
applying the NHPM to Eq.(19), we get 
1 
wn |*°Ds0(u,2) = -5(@), + 0-0], (3) 


by using the inverse natural transform to both sides of Eq.(33) and the initial condition, 


(34) 


uns 1 
O(u,t) = p2+N7 (: ree s(-) )x[-3@%, ~@2 + o| 
by applying homotopy permutation method on Eq.(34), 


S 
Dio Pry = We — N+ (2 5 +5(2) ) MER o" An — EmoP"Ba + Leao "rl, (35) 
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by comparing both sides of the Eq.(35), the following result is obtained, 
0 a = 
p :P,=e, 


S 
p:0,=N7! (2 -s+s(*) ) Mp Ay ~ pBy + p*sI} 


MAY 
p2:®, = N7 (: -~S+S (—) arto ~ p'B, + p'®,] 


By the above algorithms, 


Do = eH 


es (1 ae aa) 


1 125 
— ec = 2 _ 952 2 
®,=e (c 25 + 5°) + (28 20 Ts4 75 aa) (36) 
and so on. 
Therefore, the series solution ®(u, T) of Eq. (28) is given by 


S T2s 


®(u,t) = en! cs =39 4.52) + (GS-28?) =o ca | (37) 
If we put S > 1 in Eq.(37) , we obtain 

2 
ors 


T 
=e — 
O(u,tT) =e 41 + rT + 1 


+) a0 SU exp-u +0. (38) 
k=0 
Ultimately, the exact solution of Eq.(28). 

O(u,t) =e #tt = (39) 

Example 3. Suppose that the nonlinear system with the FABCO: 
D?@(u,t) — Pp P+ © =0,0<5 <1, 


DIW(u,t)-O, +¥+2=0,0<y<1, (40) 
where 0 < S,4 < 1 with the initial condition 

®(u, 0) = sinh(u), 

W(u, 0) = cosh(). (41) 

Below we present the NVIM, 

We obtain using the method's algorithm 


U S 
@ni1(u,T) = O(u,0) — N71 ( =545 (-) ) 0m ae a) 


Wo i(,T) = P(u, 0) + N72 ((: eo (“)) N{Oun — Pn - ®,}), (42) 


Now we discover the approximate solutions as follows: 


®, = sinh(u), ¥> = cosh(w), 


© 2023,CAJMTCS | CENTRAL ASIAN STUDIES www.centralasianstudies.org ISSN: 2660-5309 


36 


CENTRAL ASIAN JOURNAL OF MATHEMATICAL THEORY AND COMPUTER SCIENCES 


S 
®, = sinh(w) — cosh() (: —S+ saacp) 


a 
W, = cosh(u) — sinh(u) (: —A+ iar) 


Th 
d= S045) 
©, = sinh(y) + ‘ ee 1) | (sinh) — cosh()) 
ro) r(S+1)° a 1) 
5 
-S(1- 8) +S0-S) 


+ “s 725 cosh(Z), 


a (s+ a “Fs +1) 
1 
GQ=9)0 = +50-) ——— 
WY, = cosh(u) + es, — 1) (cosh(u) — sinh(w)) 
TA) Tae *"TE FATD 
—A(1-“* +41 - Dap 


+ i a sinh(). 


12 ——__ ——— 
rAt+ 7rd TGR + 1) 
Hence, Eq.(40) has approximate solution is given by 


A 
(1-s\1-D +105) 


© = sinh(y) + , pa 1) | (sinh) — cosh()) 
Peis “resis D 
S 
si=s)45055)— —— 
+ as oN 1) cosh() — --:, 
So ee 
S41) °° Tas) 
5 
(14-S\—-H +501-)—— 
Y = cosh(u) + Z ee 1) (cosh(w) — sinh(w)) 
ACS?) Tea nS 4a 41) 
i 
-x(1 —2) +40 -D —— 
+ . YD) sin(u) — +. (43) 


—72_ 2 > 
MTG+D ** FaKeD 


If we put S — 1 and A— 1 in Eq.(43), we recreate the problem solution as below. 


O(u,T) = sinh(p) (1 + - + ) — cosh(p) (x + £ +e , 
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W(u,t) = cosh(p) (t+£+--)—sinh(y) (1+54+--). (44) 


Now, we present the NHPM, 
applying the NHPM to Eq.(40), we get 


Feo wM@y(u,2) = sinh(u) — pnw [(1-5 +5 (2) ) WEY —¥ 


Fido p"Wa (ut) = cosh(u) — pw-*[(1-244(2) ) Mun 


on comparing both sides of the (45) , 
p?: 9, = 1 (x,0), p®: By = V(x, 0), 


By the above algorithms, 
®, = sinh(u), ¥5 = cosh(w), 


®, = —cosh(u) (1 —S+S a): 


WY, = -sinh() (1- A+ is), 


Th 


(1-S)(1- +20-5) aay 


®, = : wax | (Sinh() — cosh()) 


Cee | ee ry ec 
+S — Mista ArStaxaD 
S 


T 
+|a-sP 4250-9 as" 
5 


> 


W, = P “sax | (cosh(x) — sinh(q)) 


TAC — 5) Fara t ARG SKE D 
TA 


+|a- a)? + 24(1 -2I) + #? 


>TO wer +1) 


Therefore, the approximate solution of Eq.(40) is given by 
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S 
® = sinh(u) — cosh(u) (: oat Sierp) 
th 
(1-S) -4) +40 -S) ——~ 
: : Ee Dl (ania) =e0shGa) 


ct 
rots) I(S + 1) + Ms aRT 1) 
Ss 


28 
+ ja — §)* + 2S(1-S) a + §? T@s+1) cosh(u) +°*, 
a 
Y = cosh(y) — sinh(y) (: —A+ maa) 
8 
d=S)0-=1) +s =) —— 
: a re) | (cosh(u) — sinh()) 


C=) Tae “TG +4 pl 


+ ja —%)* + 2A(1 —- 4) ——~ sinh(w) + «+. (46) 


a 

+ §2 ———— 

rat Dt Ten a — 
If we put S — 1 and A— | in Eq.(46), we get 


®(u,T) = sinh(p) (1+ +5 + ) — cosh(y) (x + - fou ), 


(4,1) = cosh(y) (1+£+--) - sinh (t+ 54--).€47) 


This solution is equivalent to the exact solution in closed form: 
®(u,T) = sinh(u) cosh(t) — cosh(u) sinh(t) , 

Y(u,T) = cosh(y) cosh(y) — sinh(p) sinh(p). (48) 
Conclusions 


In this work, two nonlinear equations and a linear system of partial differential equations were successfully 
approximatively solved utilizing the NVIM and NHPM employing the FABO sense. The examples show that 
the NVIM and NHPM conclusions are very similar. The approaches are particularly successful in 
analytically and numerically solving a wide range of classes of linear and non-linear fractional differential 
equations. 
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